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We establish the jump conditions for the wavefunction and its derivatives through the 
formal solutions of the wave equation. These conditions respond to the requirement of 
continuity of the perturbations at the position of the particle and they are given for 
any mode at first order. Using these jump conditions, we then propose a new method 
for computing the radiated waveform without direct integration of the source term. We 
consider this approach potentially applicable to generic orbits. 
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1. The jump conditions 

The complexity in assessing the continuity of the perturbations at the position of 
the particle in the Regge- Wheeler gauge has led, among other motivations, to work 
in the Lorenz gauge, at the price of loosing the availability of the wave equation. 
Nevertheless, it has been indicated by two different heuristic arguments^'^ that the 
even metric perturbations for radial fall should belong to the continuity class 
at the position of the particle, in the Regge- Wheeler gauge''. Herein, we require 
that the perturbations are C'^ by identifying the conditions that the wavefunction 
and its derivatives have to satisfy for allowing the perturbations to belong to such 
category. Our analysis is based on the solutions of the Zerilli equation, not on the 
equation themselves. The inverse relations for the perturbation functions K, H2, 
Hi are given by (having suppressed the I index and being m = 0): 

K= — — ^- 1 M'r- ^ 



r2(Ar + 3M) V ^7 (A -M)(Ar 3M)r 



"The discontinuities of the wave function and its derivatives were also addressed elsewhere.^ 
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9M3 + 9AM2r+3A2Mr2 + A2(A+l)r3 SM^-AMr+Ar^ 
= rHXr+3Mr ^+ KAr+3M) ^ ^r+{r-2M)^ 

KU°{r-2M){X^r^ +2XMr-3Mr + 3M^) KU'^{r-2Mf , 
r(A + l)(Ar + 3M)2 ~ (A + l)(Ar + 3M) 

Ar^ - SMAr - 3M ^ k ^7° z„ (Ar + M) k U° z„ r(r ~ 2M) 

^ ~ (r - 2M) (Ar + 3M) (A + l)(Ar + 3M) ^ (A + l)(Ar + 3M) 

where A = l/2(/ - 1)(; + 2), k = 'im^/(2l + T)^^, 6 = S[r~Zuit)] and 5' = 5' [r-z„(t)]; 
U° = E/{l-2M/zu) is the time component of the 4-velocity, z„ the coordinate time 
dependent position, E the energy of the particle. Since the wavefunction ^ belongs 
to the continuity class, it and its derivatives can be written as, 

r) = ^+{t, r) Oi + r) 62 *,r = *+ei + ^^Oz + (*+ - S 

^,rr = *+,ei + *;,e2 + (*+ - s+{^+- f -) 5' 



^,tr = *^.ei + *;^,e2 + (*+ - f;^) ^ - (*+ - |r=.„ i„<5' 

where ©1 = © [r — 2;u(f)], and ©2 = © [zu{t) — r] are two Heaviside step distribu- 
tions. For the second derivatives, the property of the Dirac delta distribution, at the 
position of the particle: f {r)5' [r - Zu{t)] = f{zu{t))S'[r-Zu{t)]-f'{zu{t))S[r-Zu{t)], 
has been used. The discontinuities of ^ and its derivatives must be such that they 
cancel whcin combined in K , H2 and Hi at the position of the particle. Aftcir re- 
placing and its derivatives in the perturbations, continuity requires that the 
coefficients of ©i must be equal to those of ©2, while the coefficients of 6 and S' 
must vanish separately. Finally, the jump conditions for ^ and its derivatives are 
found (the jump conditions provided by K, H2 and Hi are equivalent): 

= ^-g^n +_ ^ [6M^ + 3MAz„ + A(A + 1)^] 

(A + 1)(3M + A0„) (A + 1)(2M-^„)(3M + A0„)2 

_^ _ _ KE[3M3(5A-3) + 6Af2A(A-3)z„ + 3MA2(A-l)4-2A2(A + l).z2] 

*,rr-*,rr - (A + 1)(2M - Z„)2(3M + A^„)3 

..,+ kEzuZu + _ kE (3M2 + 3MXzu - Xzl) i„ 

•* '* {2M - z^){3M + Xzu) (2M-0„)2(3M + A0„)2 
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Fig. 1. One of the possible four cases that a cell is crossed by the particle. 

2. The new integration method 

The key idea of this method is to use the above jump conditions (rewritten in terms 
of and t) rather than integrating the source term over the cells which are crossed 
by the world line of the particle. Instead, for the cells never crossed by the world line, 
the integration method retains the classic procedure. ^'^ We define a — 'i>{t + h,r^,), 
P = ^'(t, — h), 7 = 'i'{t — h,r^,) and S — '^{t,r^, + h) as the four vertices of the 
diamond of Fig. 1 centered on a = tp{t,r^,). For the case of Fig. 1, the world line 
crosses the line joining /3 and 5 at the point a and the line joining 7 and a at the 
point b. We get a set of 8 equations relating ^'+, , *+^, 

a, /3, 7, S, a, ea, £&, h, from which we deduce the value of a: 

There are other three different ways that the particle may cross the cell and 
similar relations may be drawn. The method is applicable to generic orbits. 
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